Primitive divisors of Lucas and Lehmer sequences, II 
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1. Introduction 



Let a and /3 be algebraic numbers such that a + (3 and a/3 are relatively 
prime non-zero rational integers and a/P is not a root of unity. The sequence 



(N 



defined by m„ = — — for n > is called a Lucas sequence. 



^. v-vn=u . 

If, instead of supposing that a + /3 G Z, we only suppose that (a + is 
a non-zero rational integer, still relatively prime to a/3, then we define the 
Lehmer sequence {un)'^^Q associated to a and (3 by 



r> . I — II n IS odd 

a — p 



a^-/3 



II n IS even. 



We say that a prime number p is a primitive divisor of a Lucas number 
Un if p divides u„ but does not divide (a — /3)^U2 . . .Un~i. Similarly, p is 
a primitive divisor of a Lehmer number Un if p divides Un but not (a^ — 

/3'^)'^U3. . .Un-l- 

Stewart [m p. 80] showed, as a consequence of his Theorem 1, that if n > 
C then Un has a primitive divisor, where C = e^^^2®'' for Lucas sequences 
and C = e^^^4®^ for Lehmer sequences. In Theorem [21 we shall obtain an 
improvement over Theorem 1 of [14] as well as decreasing the size of C. 
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In an earlier article [16], we enumerated all Lucas and Lehmer sequences 
whose n-th element has no primitive divisor for certain n < 12 and all 
12 < n < 30. We also presented some evidence to support the conjecture 
made there that for n > 30, the n-th element of any Lucas or Lehmer 
sequence always has a primitive divisor. 

Here, we present some further results concerning this conjecture. Our 
main result. Theorem [H states that the conjecture is true if the absolute 
logarithmic height of a is small. In addition to providing further evidence for 
the validity of the conjecture (or at least not providing a counterexample), 
this result will also be useful in a forthcoming work where we shall make 
further improvements to the size of C. 

Throughout this paper, we shall denote by h(a), the absolute logarithmic 
height of the algebraic number a, which we shall define via its relationship 
to the minimal polynomial of a over Z. Suppose that 

d 

adX" + ad-iX''-^ + . . . + ao = arf J] (X - a,) G Z[X] 

1=1 

is the minimal polynomial of a over Z with > 0, then we define 

^ logOd + X^f^ilogmax(l, \ai\) 

Theorem 1. Suppose a and /3 generate a Lucas or Lehmer sequence with 
h{f3/a) < 4. Then, for all n > 30, the n-th element of this sequence has a 
primitive divisor. 

We prove this result by using Stewart's idea [El Section 5] of looking at 
certain Thue equations. For any Lucas or Lehmer sequence (Mn)^o' there 
is a pair of integers (p, g), dependent only on the sequence, such that if Un 
has no primitive divisor then {p, q) is a solution of one of certain finitely 
many Thue equations associated to n. We use this to show that if, for 
n > 30, is without a primitive divisor then n must be the denominator 
of a convergent in the continued-fraction expansion of arccos(p/(2g))/(27r). 
The advantage gained by this is that the convergents of real numbers grow 
quite quickly and so the problem of checking each n less than 2 ■ 10^° is 
reduced to checking no more than fifty such n. 

In fact, we will show that if m„ has no primitive divisor, then the conver- 
gent k/n must be an extremely good approximation to arccos(p/ (2g))/ {2tt), 
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SO good that except for a few exceptional cases with n small, we can show 
directly that k/n is not sufficiently close to the number in question and 
therefore, eliminate n from consideration. In the remaining cases, a direct 
examination of m„ proves our desired result. 

Stewart's upper bound for n, stated above, is quite large and would thus 
give rise to extremely long calculations just to determine the convergents. 
Fortunately it is now possible to reduce this upper bound considerably. Be- 
cause of its benefit to our work here, we shall determine such a smaller upper 
bound. In fact, we establish a more general result which is an improvement 
over Theorem 1 in Stewart's paper p3], whose proof requires little more 
effort than proving the more specific result which only applies to Lucas and 
Lehmer sequences. 

Theorem 2. (i) Suppose a and /3 are algebraic integers with P/a having 
degree di over Q, (a, /3) = (1) and P/a not a root of unity. Then there is a 
prime ideal V which divides the ideal (a" — /?") but does not divide the ideals 
fori <m < n for all n > max{2(2'^i - 1), 4000(rfi log(3c/i))^2 j _ 

(ii) // a and /3 generate a Lucas or Lehmer sequence then the n-th element 
of this sequence has a primitive divisor for all n > 2 ■ 10^^ . 

2. Preliminary Lemmas to Theorem [2] 

We shall first require a lower bound for linear forms in two logarithms. 
The work of Laurent, Mignotte and Nesterenko |B] will be suitable for our 
needs. We also need a good lower bound for the height of a non-zero alge- 
braic number which is not a root of unity. 

Lemma 1. Suppose that is a non-zero algebraic number of degree D > 2 
over Q which is not a root of unity. Then 

"<-'> > fl(log(3Z)))3- 
Proof. This is Corollary 1 of □ 

Now let us continue. 

Lemma 2. Let be a non-zero algebraic number of degree D over Q which 
is not a root of unity and let log 7 denote the principal value of its logarithm. 
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Put 

A = bi log(-l) - 62 log 7 = biTci - &2log7, 
withbi a positive integer, 62 a non-negative integer and B = max(|6i|, I62I52). 
//A ^ then 

|A| > exp (-81.9(Dlog(3D))=^/i(7)(log5)2) . 

Proof. First let us suppose that I7I 7^ 1. We can write 7 = re*^ where r > 
and — 7r < 6 < 71. Since = 7 ■ 7, we have 

2h(r) = h(r2) = h(7 ■ 7) < h(7) + h(7) = 2h(7), 

so h(r) < h(7). Thus, by Liouville's inequality we have 

|A| = Ib^in - b2i9 - &2logr| > |logr| > 2"-^ exp (-Dh(7)) , 

and the lemma follows by Lemma 1 and the fact that h(7) > log 2 if = 1. 

We now turn to the case of I7I = 1. Since 7 is not a root of unity, D > 2. 

To obtain our lower bound for |A| in this case we will use Theoreme 3 
of |8]. However, this result requires that 61 and 62 be non-zero, so we must 
deal specially with the case of 61 =0. 

By Liouville's inequality and Lemma [H 

|A| > |62log7| > |log7| > 2"^exp(-Dh(7)) > exp {-2D^(log{3D)Yh{-f)) . 
It is now clear that the lemma holds in this case. 

To obtain a good constant in our lower bound we show that we may 
assume B > 679000. From Liouville's inequality, we obtain 



|A|>2-/^exp(-^?^|M) 



We can use D/2 here instead of D since 7 ^ M (see Exercise 3.4 of [18j). 
So the lemma is true whenever 

81.9Z^^(log(3D)f(log5r-|>i^. 
Since D > 2, applying Lemma [H this inequality holds if 

0.02174 

B B 

Using Maple, one can check that this is true for 2 < B < 679000. 
We now invoke Theoreme 3 of [8]. Let a = max {20, 12.85| log7| + Dh.{'~f)/2} 
and H = max {17, D log(6i/(2a) + 62/(25.77r))/2 + 2.3D + 3.25}. Then 

(1) log|A| > -9ai7^ 
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Since a > 20 and l/(2a)+l/(25.77r) < 0.0374, H < max {17, (D/2) log5 + 0.657L> + 3.25}. 
Moreover, B > 679000 implies that {D/2) log 5+0.657^+3.25 < 0.66994L)log5. 
As this last quantity is greater than 17 for B > 679000, we have H < 
0.66994L)log5. 

We also want an upper bound for a in terms of D and h(7). First notice 
that I log7| < 71. Therefore, 12.85| log7|+Dh(7)/2 < Dh{-f){A0.37/{Dh{-f)) + 
1/2). Since D > 2, we can apply Lemma [TJ We obtain 40.37/(Z^h(7)) + 
1/2 < 20.185(log(3L'))3 + l/2 < 20.272(log(3L'))3. Therefore, 12.85| log7| + 
L'h(7)/2 < 20.272(log(3L'))3Dh(7) for all D > 2. Moreover, this quantity 
is greater than 20, so a < 20.272(log(3D))3Dh(7). 

Applying these estimates to ([1]), we find that our lemma holds. □ 



Suppose that a and (3 are algebraic integers in a number field IK of 
degree d over Q. Letting Ki = Q(/3/a), a number field of degree di over 
Q, we set (3 /a = /3i/ai, where ai and f3i are algebraic integers in Ki and 
= Ai. We may assume, without loss of generality, that |ai| > \f3i\. 

We note that, unless we state otherwise, log 2; shall always denote the 
principal branch of the logarithmic function. 

Now let us prove: 

Lemma 3. (i) We have 

log2 + log lail > log \ai - > log |q;i| - di{h{/3i/ai) + log2). 
(ii) For di >2 and n >2, we have 

log2+nlog|ai| > logK > nlog |ai|-81.97(rfi log(3rfi))^/i(/3i/ai)(logn)l 

Proof, (i) We can write ai — /3i = a;i(l — /3i/ai). By Liouville's inequality, 
log |/3i/«i - 1| > -rfi(log2 + h(/3i/«i)), 

and the result follows. 

(ii) The upper bound follows directly from the triangle inequality and 
our assumption that |q;i| > 

For the lower bound we write the quantity in question as 

nloglail +log|(/3i/ai)" - 1| . 



Applying Lemma 2.3 of with r = 1/3 and z = nlog(/3i/ai), we see 
that either 

- 1| > i or |A| = |nlog(/3i/ai) - 2k7Ti\ < 1.3 - 1| < 0.5. 

3 

In the first case, the lemma holds so we need only consider the second 
case. Here, we must have 

|Im(rzlog(/3i/ai)) - 2fc7r| < 0.5. 

Since we took the principal value of the logarithm of we have 

-TT < lm(log(/3i/ai)) < Tc and so \k\ < n/2 + 0.5/(27r) or \2k\ < n. 

As f3i/ai is, by assumption, not a root of unity, A 7^ and, since n >2, 
we may apply Lemma [2] giving 

|A| > exp (-81.9(rfilog(3di))^h(/3i/ai)(logn)2) . 

By Lemma dl we have log (1.3) < 0.07(((ii log(3(ii))^h(/3i/a;i) log^ n, since 
di >2 and n >2. Our lemma follows. □ 

Lemma 4. Let $„(X,r) = y^^^VnW^) where </)„(X) is the n-th order 
cyclotomic polynomial. Suppose that V is a prime ideal in K which divides 
for n > 2(2'^^ - 1). This implies that V divides (a" - /?"-). //, in 
addition, V divides (a™ — for some m < n, then 

ord73$„(a, /3) < ord-pn. 
Proof. This is Lemma 4 of [13]. □ 

Finally we need to bound some arithmetic functions which will appear 
throughout this article. 

Lemma 5. (i) Let uj{n) denote the number of distinct prime factors of n. 
For n > 3, 

1.3841 log n 



log log n 

(ii) For n > 3, 

n 



e'l log log n + 2.50637/ log logn' 
where 7 = 0.57721 . . . is Euler's constant. 

Proof, (i) This follows from Theoreme 11 of [TT]. 

fii) This is Theorem 15 of 112]. □ 
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3. Proof of Theorem [2 

We may also assume that di > 2, for otherwise we can write a = /3ci/c2 
where Ci, C2 G Z with (ci, C2) = 1 and so a" — /3" = (c" — C2)(/3/c2)". Now 
Zsigmondy [20] and, independently of him, Birkhoff and Vandiver [2J have 
shown that for n > 6 the n-th element of such sequences always has a 
primitive divisor. 

Therefore, we may also assume that n > 3900(2 log(3 • 2))^^ > 1.74 ■ lO^'', 
since Theorem 2 does not apply for smaller n when di > 2. 
We note that 

<l>„(a,/3) = = /3'^W<l>„(ai//3i, 1) = {(3 / (3,r^^^<!>^{a,, (3,). 

Letting A be the extension of Ai in K, we have = A^^, since 

(a, /3) = (1), and so 
(2) 

(di/d) log \Nk/q ($„(a, = log \Nk,/q ($„(«!, | -</^H log A^^./q (^1) . 
Since 

$n(ai,/5i) = n(«r-/9rr^"/'"\ 

m|n 

the right-hand side of ([2]) is 

where Moo{Ki) denotes the set of all archimedean absolute values defined 
on Ki up to equivalence. 

Applying Lemma [31 we see that the inner sum in the first term of this 
expression is at least 

\og{max{\ai\,„\(3i\y)}^n{n/m)m- log2 

m\n m\n,m>l 

li{n/m)=—l 

-81.97(rfi log(3rfi))=^h(/3i/aO ^ (logm)^ - di(h(/3i/aO + log2). 

m I n , m > 1 

^(n/m)=l 

Combining this lower bound with 

logmax(|ai|^„ \(3i\y) - logA^^i/Q {Ai) = h{(3i/ai) 
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and 

^^m/i(n/m) = (p{n), 

m\n 

we obtain 

(di/ci)log|iVA7Q($„(a,/3))| > ^(n)h(/3i/ai) 

-81.97rft(log(3di))=^h(/3i/ai) ^ (logm)^ 

m I n 

/j,(n/m)=l 

- 5^ rfilog2-rf2(h(/3i/ai)+log2). 

m I n 

^(n/m)=— 1 

Notice that n has 2'^*^"^""^ factors m which satisfy ^(n/m) = 1 and the 
same number of factors m satisfying n{n/m) = —1. Now, by Lemma [Hand 
our lower bound for n, 

d?(h(/3i/ai)+log 2)+ J2 ^1 log 2 < 0.005-2^(")dt(log(3rfi))=^h(/3i/ai) log^ n. 

fj,{n/Tn)=—l 

Thus 

(4) 

(rfi/rf) log |iV;,/Q >V;Hh(/3i/«i)-2'^(")40.9Mt(log(3rfi))3h(/3i^ 

for rfi > 2 and n > 1.74 ■ 10^°. 

By LemmaHl if \Nk/q > 't-'^, then there exists a prime ideal V 

which divides (a" — but does not divide (a™ — /?"*) for any m < n. Using 
(jlj) and Lemmadl as well as our assumptions that di >2 and n > 1.74-10^°, 
this condition is satisfied if 

> 41dt(l„g(34))=. 

From Lemma [5l we find that 

V^H . „0.3495 

2'^W(logn)2 
for such n. Therefore, ([5]) is satisfied for 

(6) n > 41200rf}^-^^(log(3t;i))^-^^ 

Since di >2, part (i) of the theorem holds. 

(ii) Let (Mn)^o ^ Lucas or Lehmer sequence generated by a and 
/3. Since a(3 and (a + are relatively prime non-zero rational integers. 
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there exist two integers p and q such that a and (3 are the two roots of 
— yjp + 2qX + q. Therefore, a, /3 = (a/p + 2g± a/p — 2g)/2 and so either 
or /3/q; is equal to (p+ a/p^ — 4g2) / (2g). Therefore we can take cii = 2 
and so part (i) of theorem imphes part (ii). 

4. Preliminary Lemmas to Theorem [1] 

Lemma 6. Let a he a non-negative real number. If x,y G M with —1 < 
x,y < 1 and \x — y\ < a then 

[a 

|arccosx — arccosyl < vrW-. 

Proof. This result follows from finding the minimum value of the function 

cos a; — cosy 

on the area in defined by < x, y < vr, x 7^ y which is 2/7r^ and then 
applying the contrapositive. □ 

Let us collect here various notations which we shall use throughout the 
remainder of this article. 

Notations. Given a complex- valued function / defined on C, we use 
to denote max|a;|=i 

For a positive integer n, we let gn{x) G Z[x] be the minimal polynomial of 
2cos(27r/n) over Z; its degree is ip{n)/2 if n > 3. We shall put Gn{X,Y) = 

F^(")/2^„(x/r). 

We let m be the greatest odd square-free divisor of n. For such m, we 
shall write hmiX) = (X™ - 

Finally, for n > 1, we let P(n) denote the largest prime divisor of n. 

As we shall see in Section 5, the crucial result needed in the proof of 
Theorem [1] is a good lower bound for |5f,^(2 cos(27rj'/n))| for {j,n) = 1. 

We will show that we need to obtain an upper bound for the absolute 
value of hm{X) on the unit circle which we find using an idea and a result 
of Bateman, Pomerance and Vaughan pQ . 

Let us start linking these two polynomials now. 
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Lemma 7. Let n>3, l<j<n with {j,n) = 1 and Cn = exp(27ri/n). 
Then 



|^;(2cos(27rj/n))| 



2sm(27rj/n) 



Proof. We can write 

MX) = n (^-a)(^-c^)= n (^'-(a+CT^ + i) 



l<i<n/2 

{j,n)=l 



l<j<n/2 

{j,n)=l 



- n + 1 - 2 cos(27rj7n)X) = ^„ (^^) 

l<7<n/2 ^ ' 



l<j<n/ 

(i,n)=l 



If y = (X^ + 1) /X then X = (y ± VY^ - 4) /2 = /(y) and so 
^ _ MfjY)) 

^ 2/(y)0U/(y))r(y) - Uf{Y))^{n)r{Y) 

yn\ ) 2/(y)¥'(")/2+i 

Since /(2 cos(27rj/n)) = cos(27rj/n) ± i sin(27rj7n), 

0^(cos(27rj/n) ± i sin(27rj/n))/'(2 cos(27rj/n)) 



^;(2cos(27rj/n)) = 



(cos(27rj/n) ± i sin(27rj7n)) 



¥>(n)/2 



Notice that /'(y) = (1 ± y/ V^"^ - 4) /2 so that 



/'(2cos(27rj7n)) 

Hence, 

|5;(2cos(27rj7n))| = 
from which the lemma follows. 



cos{2'n-j/n) 
sin(27rj/n) 



□ 



To work with the cyclotomic polynomials we shall need some relation- 
ships which they satisfy. We give these in the next lemma. 

Lemma 8. (i) Let n be a positive integer and let m be its greatest odd 
square-free divisor. We put m' = gcd{2,n)m. Then 
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(ii) Let p be a prime number and n any positive integer not divisible by p. 
Then 

4>n(X) 

(iii) Let m,m' and n be as above. We put n' = n/ gcd{n,2), hm{X) = 
^X"^ — 1)/0„(X) and Cn = exp{2TTi/n). Then, for all j with {j,n) = 1, we 
have 



n' 



Proof, (i) This assertion follows easily from the two relations: 
MX) = M-X) and 0„(X) = (^X"/™') , 

which are parts (iv) and (vi) of Proposition 5.16 from Chapter 2 of Karpilovsky's 
book [7]. 

(ii) This is again from Proposition 5.16 from Chapter 2 of [7]. 

(iii) Applying part (i), we find that 

" " m' 
NowX--l = /i„(X)0„(X) somX— 1 = hUX)MX)+h'^{X)MX). 
Letting X = (— l)™'^^Cm'' which is always a primitive m-th root of unity, 

we have (-1)(— DK+DmC^""'^ = /im((-l)"'+^C.)0™((-l)'"'+'C') and 
the result follows. □ 

We see now that we have reduced the problem of bounding \g'^\ from 
below for primitive n-th roots of unity to bounding \hm\i from above. To 
deal with this new problem, we shall now use ideas from [T]. 

Lemma 9. Let m = pi . . .p^ where pi,P2, ■ ■ ■ ,Pk o,i"g odd primes arranged 
in increasing order. Then 

k-i 

\hUX)\, < 2l[pf-'-\ 

i=l 

In fact, if k > 3 then the factor of 2 is not needed. 
Proof. From Lemma [8]^ii) , 

(7) hux) = ~ ^^^'rt'J^^ = K....v.-AX'^)M...^AX). 

Vpi...Pk~i\^ ) 

We now use induction on k to prove the lemma. 
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Since hi{X) = 1, V(X) = X - 1 and hp^p^{X) = (X^^ - l)0pi(X), the 
lemma is true ioi k < 2. 

For k = 3, we have j/ipipapsli — I^pip2lil0pip2li- Using the result just 
established for k = 2 and a theorem of Carlitz [3] which shows that |0pip2 li < 
P1P2/2, we obtain j/ipipapali ^ PiP2- This is the desired inequality for k = 3. 

Suppose now that the lemma holds for some A; > 3. We apply the 
following estimate of Bateman, Pomerance and Vaughan, which follows from 
Theorem 1 of their paper [1] and holds for A; > 3, 

k-l 
i=l 

Thus from ([7j), we have 

k-l k-l k 

\Ki---Pk+i\i < \Ki---Pk\i\(f>Pi...p,\i < [[Pi ^Pk[[Pi = [[Pi ■ 

i=l i=l i=l 

Hence the lemma holds. □ 

We need the next lemma to deal with the case q = 2, although we will 
use it for all g. A simple application of the triangle inequality would quickly 
yield the inequality below with 3|g|/5 replaced by \q\/2. However, in the case 
of g = 2, this would not be sufficient to prove our theorem: with the lower 
bound that the previous lemmas imply for \g'^{2 cos(27i k/n))\, the upper 
bound we would obtain for the left-hand side of ffTTl) would not decrease 
with n but actually grow with n. To refine this trivial estimate, it seems we 
must resort to an argument like the one which follows. 

Lemma 10. Let n > 30 be a positive integer and let p and q be non-zero 
integers with q > 2, \p\ < 2q and 

\Gr,{p,q)\<P{n/{n,3)). 

For 1 < j < n/2 with (j, n) = 1, we put I3n^ = p — 2qcos{27rj/n). Define 
k by \f3i^^\ = min \f3i^^\. Then 

j^\...n/2 

{hn)=l 

(3|g|/5)^(")/2|^;(2cos(27rfc/n))| 



n A 

j¥'k,U,n)=l 
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Proof. We first divide the interval (—2, 2) into four subintervals and divide 
tlie set of integers less than n/2 which are relatively prime to n into four 
associated subsets. Let A = (—2, —1),^' = {m : n/3 < m < n/2, {m,n) = 
1},B = (-1,0), -B' = {m : n/A < m < n/3, (m,n) = 1},C = (0,1), C = 
{m : n/6 < m < n/4,{m,n) = = (1,2) and V = {m : < 

m < n/6, (m, n) = 1}. If we let !f{k,q,n) denote the number of inte- 
gers in the interval {nq/k,n{q + l)/k) which are relatively prime to n then 
l^'l = ip{6, 2, n) = ip{n)/2 - ^{3, 0, n), \B'\ = ip{3, 0, n) - ip{4, 0, n), \C'\ = 
(/9(4, 0, n) — ip{6, 0, n) and I'D'] = ip{6, 0, n). 

Using Theorems 5-7 of [9], we have the following inequalities for the 
cardinalities of these sets of integers: 



ip{n) - 2^(") 



(8) 





6 






-3- 


2w(n) 




12 




(f{n) 


-4- 






12 




(f{n) 


- 2 • 
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< l^'l < 

< \B'\ < 

< \C'\ < 





+ 2^ 


(n) 




6 




(f{n) 


+ 3 






12 




(p{n) 


+ 4 






12 




(p{n) 


+ 2 
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Let us observe that p/qEAUBUCUV and 



(9) 



n/2 

n 



jytk,{j,n)=l 



n/2 

n 



1 - 



-1 



j^k,{j,n)=l 



|^?;(2cos(27rA;/n))|. 



If p/q G >4. then p < —3, since q > 2, and either /c G or A; is the largest 
element of B'. Thus, /3n ■* > 3 + 4cos(27rj/n) for each j G C U and so 



I <ci= (P(n/(n,3)) J] (3 + 4 cos(27rj7n))-M 
\ j&c'uv / 



l/(An)/2-\C'\-\V'\) 



Combining these inequalities with ([9]), we obtain 



n/2 

n 



j¥'k,{j,n)=l 



I |^(n)/2-l 

n 



2l^'| + |B'| 



1 + 



Cl 



3 + 4cos(27rj/n) 



|(7;(2cos(27rA;/n))| 
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Now suppose that p/q G B. If /3f ^ < then ^^^^ < for j G C U V, so 
|1 - I3n^/I3n \ < 1 for such J and 



n/2 

n 



jjtk,{j,n)=l 



\ Mn)/2-l 

^ 2\A'\+m l^;(2cos(27r/c/n))|. 



Since the quantity before |5'^(2 cos(27rA;/n))| on the right-hand side of this 
expression is at least as large as the similar quantity obtained ior p/q e A, 
we can ignore this case. 

If plr^ > then |1 - pi^^/pi^^l < 1 for j e A' so a similar analysis to 
that above shows that 

^ l/{v(n)/2-\V'\) 

I <C2= (P(n/(n,3)) J] (1 + 4 cos(27rj7n))-^ ' 



and 



n/2 

n 



i=i 
j7^k,{j,n)=l 



> 



2\B'\+\C'\ 11 I 
jev 



C2 



1 + 4 cos(27rj7n) 



|5;(2cos(27rfc/n))|. 



11 p/q e C then, by the same reasoning, we obtain 

< C3 = ( P(n/(n,3)) J] (1 - 4cos(27rj7n))-^ 



l/{^{n)/2-\A'\) 



and 



n/2 

n /^i^'^ 



i=i 
jVfe,(j,n)=l 



> 



|g|.(n)/2-l / 

2|B'HC'I llA l-4cos(27rj7n) 



-1 



If p/g e X>, then 

I <C4= [P(n/(n,3)) n (3-4cos(27rj7n))-^ 



|5;(2cos(27r/c/n))|. 



l/(¥>(n)/2-|.A'|-|B'|) 



and 



n/2 

n 



j^k,{j,n)=l 



2|C'|+p'| 



C4 



3 — 4cos(27rj7?T') 



|^;(2cos(27r/c/n))|. 
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For n < 210, n = 231 and n = 462, we can use these estimates to show 
by direct calculation that our lemma holds. 

To deal with n > 210, we first show that max(ci, C2, C3, C4) < 1 for such 
n. Using the above expressions for these quantities we see that this holds if 
min(3l^'l,3l^'l) > n. By dH]), both \A'\ and \V'\ are at least ¥?(n)/6-2'^("V3, 
so we need only prove that {(fin) — 2 ■ 2'^^"))(log3) > 61ogn for n > 210. 

For 210 < n < 330 = 2 ■ 3 • 5 ■ 11, 2'^(") < 8 < n°-'^^^. Lemma [5](ii) yields 
the lower bound ip{n) > n'^-'^^^ for n > 210. Since logn < n^-^^^ for n > 210, 
we need only show that 77,0 075(^^0.33 _ 2) log 3 > 6 for n in this range. But 
this is easily seen to be true. 

For 330 < n < 2310 = 2 ■ 3 • 5 ■ 7 • 11, 2'^(") < 16 < n°-^^. Moreover, 
by Lemma EKi), for n > 2310, 2'^(") < n^-^^^^f ^ogiogn ^ ^o.47_ Therefore, 
for n > 330, 2'^^"^ < n'^-^^. Applying Lemma E](ii) again, we find that 
(p{n) > 77, "'^^ for n > 330. Since logn < n^'^^ for n > 330, we need only 
show that 77,0.17^^^0.25 _ 2) log 3 > 6 for n in this range which is also easily 
seen to be true. Therefore, max(ci, C2, C3, C4) < 1. 

So, from our lower bounds for the absolute values of the products of the 
given above, to prove the lemma we need to show that 

max (2\^'\MB'\^^fr^ynMT>'\^2\^'\Mn^ 

is less than (5/3)^('^)/2_ 

Let us first show that 2\^'\{2/3)\^'\ > 1 and 2l^'l(2/3)l^'l > 1. These in- 
equalities will show that either the first or the last terms give the maximum 
in this expression. 

For the first of these two inequalities to be true, by ([HD we need to 
show that 0.08ip{n) — 0.26 ■ 2'^^'^^ > 0. Similarly, the second inequality 
requires that the stronger inequahty 0.08(/?(n) - 0.34 ■ 2'^^'') > holds. So 
we need only consider this last inequality which we shall rewrite in the form 
0.08/0.34 > 2'^("V^H- 

For 210 < n < 330, we saw in a previous paragraph that 2^^'^'^ / Lp{n) < 
^-0.33 < 0.171 < 0.08/0.34. We also saw that 2'^(")/(/?(n) < < 0.235 < 

0.08/0.34 for n > 330. Therefore, our desired inequalities holds and we need 
only try to bound 2l^'l+ie'l(4/3)l'^'l+l^'l and 2^^' 1+^(4/3)1 A' \+\b'\ f^^^ 

Notice that \A' U B'\ = (p{A, 0, n) and that \C' UV'\ = y?(4, 1, n). Lehmer 
[U p. 351] has noted that £'(4, 1, n) = —£'(4, 0, n), where E{k, q, n) denotes 
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(p{n) — k(f{k,q,n), so we need only examine 

(8 /3)'^^"^/'*(3 /2) l^(4.0:")l/4_ 

Lehmer also gives precise information about E{A,0, n) in Theorem 6 of 
[9]. If n > 4 and 4 divides n or n is divisible by a prime congruent to 1 
mod 4 then i?(4,0,n) = and our proof is complete. If neither of these 
conditions is true then |i?(4,0,n)| = 2"^*^"'^ where n' is as in the statement 
of Lemma [Hl(iii). Notice that when i?(4,0,?7.) 7^ 0, n is not congruent to 
mod 4, so n' is the odd part of n. 

A direct calculation shows that for 210 < n < 750, with the exceptions 
of n = 231 and 462 which we considered above, |i?(4, 0, n)|/(y9(n) < 0.05. 
Therefore, (8/3)^(")/^(3/2)l^(^'0'")l/^ < (5/3)^(")/2 210 < n < 750, n ^ 
231,462. Recalling that we showed by calculation that the lemma holds for 
30 < 72 < 210, for n = 231 and for n = 462, we now know that the lemma 
holds for 30 < n < 750. 

Notice that if n < 4389 = 3-7- 11-19 then either \E{A,0,n) \ = or 2'^("') < 
8, since in the latter case n' is odd and without prime divisors congruent to 
1 mod 4. Using Lemma [5t^ii) , ip{n) > 160 and so \E{4,0,n)\/(f{n) < 0.05 
for n > 750 and our lemma holds for 30 < n < 4389. 

Applying the inequality 2^^^"^ < ^0-96/ log log which follows from LemmaEI^i), 
and part (ii) of this same lemma, we find that 

2<^(n) /4 ^0.96/ log log « ( 1 . 78 1 1 log log n + 2 . 5 1 / log log n) 

The right-hand side is a monotone-decreasing function for n > 10 and 
so it is less than 0.05 for n > 4389. Therefore (8/3)^(")/^(3/2)l^(^'0'")l/^ < 
(5/3)'''^")/^ for n > 4389, which shows that the lemma is true. □ 

5. Proof of Theorem [1] 

Let (ttn)^o ^ Lucas or Lehmer sequence generated by a and (3. As 
noted in the proof of Theorem |2l there exist two integers p and q such that 
a and /3 are the two roots of X"^ — y^p + 2qX + q. Notice that the n-th 
element of the sequence generated by ia and iP is just ±m„. Therefore, we 
can assume that q = a(3 is positive. Also notice \p\ < 2q for otherwise a 
and P are real and Carmichael [1] , Ward [19] and Durst [5j have shown that 
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in this case the n-th element of these sequences has a primitive divisor for 
n > 12. 

Let us define the /3n ■''s and pit^ as in Lemma [TUl Stewart [HI Section 5] 
has shown that if the n-th element of this sequence has no primitive divisor 
then 



(10) 



\Gn{p,q)\= n W^\<Pin/{n,3))foin>U. 



l<j<n/2 

(i.")=i 



Since \p\ < 2q, upon applying Lemma [TUl we obtain 

, P(n/(3,n)) ^ (5/3)^(")/2p(n/(n,3)) 

\' n — -1— r I „ / . % 



Y[ \(3l^^\ ^ |^;(2cos(27rA;/n))||g|^(")/2-i^ 



l<j<n/2 

j;VA:,(jr»=l 

for n > 30. 

Therefore, if we can show that 

'2nk 
q \ n 

then, by Lemma [6l 



:ii^ 



2 cos 



, (5/(3|g|))^(")/^P(n/(n,3)) ^ 4 



\g'^{2cos{27rk/n))\ 



1 

2^ 



arccos 



2q 



< 



2n2' 



and so, by Theorem 184 of [6J, k/n must be a convergent in the continued- 
fraction expansion of arccos(p/(2g))/(27r). 

Hence we first want to show that for n sufficiently large, the right-hand 
inequality of flTTl) holds. We start by considering the case of g = 2, as this 
is the most difficult one. 

5.1. The case q = 2 

Using the notation of Lemmas [8] and [9], we find, from Lemma [H that 

for m > 1. 

If m = 1, but n > 1, we have \hmiX)\i = 1 < n2'^*"^"V^(«). 
Combining this upper bound with Lemmas [7] and [^iii), we obtain 



|(7;(2cos(27rA;/n))| > 



n 



4^2"(")-Va;(n) 



for n > 1. 
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Applying this lower bound to the right-hand inequality of ( ITTj) and squar- 
ing both sides, we want to show that 

for n > 30. 

To prove that this holds for n sufficiently large, we take the logarithm 
of both sides, which yields 

¥p(n)log(5/6) + 2'^(")(logn)/w(n) + 81ogn < 0. 

From Lemma 0(11), we see that f{n) > n^-^^'^^ for n > 3500, while for 
the term involving u{n) we use the fact that is a monotone increasing 
function for x > 1/ log2, 2^/1 = 2^/2 and Lemma[5](i). In this manner, our 
problem is to show that 

„0.96/loglogn Incrlnp-r) 

-0.182n™ + "1 l^Si^ + 81ogn < 0. 

1.384 ^ - 

For n > 3500, the sum of the second and third terms is at most nP-^^^'^, 
Therefore, we need only show that — 0.182n'^'^°^ + 1 < 0, but this is easily 
seen to be true for n > 3500. So we have an initial bound of intermediate 
size. 

Notice though that we did not make full use of the Lemma |9] in this 
argument. A direct calculation on a computer using the result given in 
Lemmas [TfHT iii) and |9] shows that the right-hand inequality of f|TT]) holds for 
all n > 1260 when q = 2. 

In the case of g = 2, Lucas and Lehmer sequences can result from p = 
—3,-1,1 and 3. Since Gn{p,q) is a product of terms of the form p — 
2qcos{2m/n), it is quite easy to calculate Gn{p,q), although care must be 
taken to maintain sufficient accuracy, and so we can check whether m„ has 
primitive divisors by means of fllOp . However, to check Un for each n up to 
1260 in this manner is quite time-consuming. Fortunately, one can quickly 
extract still more information from f lTT]) . Given n,p and q, it is easy to 
find the integer k with {k,n) = 1 which minimizes the far left-hand side of 
fllip . As when considering 1260 < n < 3500, we can bound from above the 
middle quantity in f lTT]) . For q = 2,p = —3, —1, 1, 3 and 330 < n < 1260, we 
can verify in this way that the left-hand inequality in ( ITT]) is violated and 
so for such n, the n-th element of these sequences has a primitive divisor. 
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But we still need to consider 30 < n < 330. For these n, we use ( ITOj) as 
described earlier in this paragraph. 

For n > 1260, we have seen that n must be the denominator of a con- 
vergent in the continued-fraction expansion for arccos(p/(2g))/(2-7r). 

The question arises of how to deal with these n. We are fortunate that 
in these cases the middle quantity in f llip is extremely small. For such n, 
we proceed in the same manner that we checked the left-hand inequality in 
(fTT]l holds for 330 < n < 1260, except that now we know k too. Theorem [2] 
tells us that we need only check those convergents k/n with n < 2 ■ 10^^. 
For each convergent computed with n < 2 ■ 10^°, \p/q — 2 cos(27r/i;/?2)| was 
considerably larger than the bound that the left-hand inequality of ffTTl) 
requires if m„ were to be without a primitive divisor. In Table 1, for p = 
—3, we list the convergents with n > 1260 and give the logarithms of the 
required and actual bounds, denoted c/req and ci^ct' respectively. The value 
of log|c?req| given in Table 1 is truncated to its integer part, whereas the 
value of log l^actl truncated to one decimal place. 

Proceeding in this same way for p = —1, 1 and 3, we are able to conclude 
that if (m„)5^o is the Lucas or Lehmer sequence generated by any of the pairs 
{a, /3) = (1 ± v^)/2, (V3 ± v^)/2, ( ^5 ± v^)/2 or (v^ ± v^)/2, 
then Un has a primitive divisor for n > 30. 

5.2. The case of g > 2 

For such pairs {p, q), we proceed along the same lines. The only difference 
is that less work is required for small n. We already know that the right- 
hand inequality of f lTT]) is satisfied for n > 1260 by our work in the previous 
section. We can check directly, as with 1260 < n < 3500 for q = 2, that the 
right-hand inequality of f llip holds for ng_i > n > Ug where Ug is given in 
Table 2. 

As in the case of 30 < n < 330 for g = 2, we directly check those Un 
with 30 < n < Ug for primitive divisors and for larger n we compare the 
required and actual differences of \p/q — 2 cos(27rfc/?7,)| in f|TT]) to estabhsh 
our result. The actual difference is less than the required difference for all 
Uq < n < 2 ■ 10^° and all 3 < g < 3000 (this corresponds to all pairs of a 
and P with h(/3/a) < 4). By Theorem [21 Theorem [T] now follows. 
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k 


n 


log dreq 


log c^act 


497 


1291 


-116. 


-12.6 


579 


1504 


-68. 


-13.7 


1655 


4299 


-260. 


-15.4 


3889 


10102 


-459. 


-18.9 


52212 


135625 


-8207. 


-22.1 


56101 


145727 


-12970. 


-22.4 


108313 


281352 


-8086. 


-24.3 


381040 


989783 


-90228. 


-26.1 


489353 


1271135 


-90181. 


-26.7 


870393 


2260918 


-93683. 


-28.3 


2230139 


5792971 


-493472. 


-29.5 


3100532 


8053889 


-734197. 


-30.8 


8431203 


21900749 


-1745895. 


-32.3 


11531735 


29954638 


-1165244. 


-33.1 


19962938 


51855387 


-3104401. 


-34.1 


31 19 1673 


81810025 


-5101005. 


-35.2 


51457611 


133665412 


-5943915. 


-35.8 


82952284 


215475437 


-19412834. 


-38.5 


798028167 


2072944345 


-144472147. 


-41.8 


1679008618 


4361364127 


-374075698. 


-42.8 


2177036785 


6131308172 


-293278281. 


-11.3 


6633082188 


17229981071 


-1438733756. 


-45.4 



Table 1. {p,q) = (-3,2) Verification 



All the calculations in this article were performed using Release 3 of 
Maple V and UBASIC 8.74 on an IBM-compatible PC with an 486DX2 
running at 66 MHz. In total, the calculations required just over 100 hours on 
this machine. Many of the calculations were performed using both systems 
to provide a check on the quantities obtained and the results were always 
identical up to the specified accuracy. 
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Q 


Ug 


2 


1260 


3 


330 


4 


210 


5 


120 


6 


90 


7 


78 


8 


66 


9,10,11 


60 


12,..., 20 


42 


> 21 


30 



Table 2. Values of Ug 
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